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ABSTRACT 



We present high angular resolution (~ 3") and sensitive 1.3 mm continuum, cyanogen (CN) and vinyl cyanide (C2H3CN) line ob- 
servations made with the Submillimeter Array (SMA) toward one of most highly obscured objects of the W51 IRS2 region, W51 
North. We find that the CN line exhibits a pronounced inverse P-Cygni profile indicating that the molecular gas is infalling inwards 
this object with a mass accretion rate between 4 and 7 x 10"" Mq yr"'. The C2H3CN traces an east-west rotating molecular envelope 
that surrounds either a single obscured (proto)star with a kinematic mass of 40 Mq or a small central cluster of B-type stars and that is 
associated with a compact high velocity bipolar outflow traced by HtO masers and SiO molecular emission. We thus confirm that the 
W5 1 North region is part of the growing list of young massive star forming regions that have been associated with infalling motions 
and with large mass accretion rates (~ 10"^ - 10"* ), strengthening the evidence for massive stars forming with very high accretion 
rates sufficient to quench the formation of an UCHII region. 

Key words, stars: pre-main sequence - ISM: jets and outflows - ISM: individual: (W51 IRS2, W51 North, G49.49-0.37) - ISM: 
Molecules, Radio Lines - ISM: Circumstellar Matter - ISM: Binary stars - ISM: Envelopes - 



1. Introduction 

One of the main questions related to star formation is whether 
massive stars (> 10 M©) are formed through gas accretion via 
a circumstellar disk/torus or whether other mechanims play a 
role. It was believed that the powerful radiation fields and stellar 
winds produced at the very beginning of nuclear burning will 
increasingly inhibit farther accretion of material thereby lim- 
iting the maximum stellar mass to about 10 Mg (lKahnlll974l: 
Larson & Starrfield 1971; Yorke & K ruegellfTOTTl) . Several the- 
oretical models have since been proposed to solve this puz- 
zle: the form ation of massiv e stars through d ense d isks with 
jets/outflows dNakanol 119891: |jinnr& Adams Il996l) . through 
merging of smaller stars ("Bonnell et al.lll998!) . through turbu- 
lent accretion (Mc Kee & Tan 2003), through competitive accre- 
tion jBon nell et al.i l2001h . and through ionized accretion flows 
dKeto & Woodi20()6h . However, due to the lack of good obser- 
vati onal guidance, these alte rnatives have remained controver- 
sial (Zi nnecker & Y orke 2007). 

With a total bolometric luminosity of about 3 x 10^ L© the 
W5 1-1RS2 region is one of the most luminous massive s tar form- 
ing regions in our Galaxy, dErickson & Tokunagal 19801). It is lo- 
cated 6-7 kpc away iri the Sagittarius spiral arm dGenzel et al.l 
119811: llmaietal .1 120021. We adopt here a distance to the W51 
region of 7 kpc. The W51 1RS2 region comprises a complex 
group of highly obscure d young objects (with n ot mid-infrared 
counte rparts, see please iKraemer et al.l d200ll) : lOkamoto et al.l 
d200l ) called "W51 North" and "W5 1d2", and a cluster 
of massive strong infrared ZAMS stars (iKraemer et al.l I2OOU 



lOkamoto et"ani2001l: iLacv et al.ir2007l) . which the most promi- 
nent member being the source "IRS2d" associated with an ex- 
tended edge-brightened cometary HII region called "W51d" 
(JGaume et al.lll993l:lLacv et alJl20(J7l) . 

The W51 North object shows strong thermal dust emission 
at (sub)millimeter wavelengths, molecular emission from a large 
(~ 4x10^ AU) hot co re at an excitation temperature of 100-200 
K dZhang et al. 19981). and very faint centimeter free-free emis - 
sion dGaume etal.lll993t IZhang et alJ[T998l: lEisner et al.ll2002l) . 
suggesting that it is forming an extremely young massive star. 
This object, in addition, contains a group of strong H2O, OH and 
SiO masers that are located in the center of this molecular and 
dusty structure (called "The Dominant Center") dSchneps et al 



Send offprint requests to: 
lzapata@mpifr-bonn.mpg.de 



Luis Zapata, e-mail: 



1981; Gaume & Mutelll 19871; iHasegawa et al.| [r986; Uki ta et al. 
1987i ; iMorita et alJll992 Y Observations of the proper motions 
of the H2O masers (which traces shocks in dust-laden gas close 
to the exciting protostars, Elitzur 1992) revealed the presence 
of a compact (~ 7000 AU) northwest- southeast high velocity 
(> 100 km s-') o utflow dSchneps et alJll98ll;lEisner et alj|2002i; 
llmai et ani2002l) . Moreover, high angular resolution observa- 
tions showed that the SiO masers s eem to be tracing t he inner- 
most parts of this powerful outflow dEisner et al.ll2002l) . Finally, 
this object has previously been identified with spectroscopic 
signatures of dynam ical collap se using emission from HCO^ 
dRudolphetal.lll99 0). and SO2 dSollinsetal.ll2004l) . 

Here we present 1.3 mm continuum, cyanogen and vinyl 
cyanide line observations toward the W51 North region made 
with the SMA. We report the presence of molecular gas accre- 
tion onto a 40 M© (proto)star or a small central cluster of B stars 
located in the center of W5 1 North region, and with an accretion 
rate between 4 and 7x 10"~ M© yr"' . 
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Fig. 1. Integrated molecular emission of the lines CNN=2-1, J = 5/2-3/2, F = 5/2 - 3/2 and N=2-l, 7 = 5/2- 3/2, F = 7/2 - 5/2 obsei-ved with the SMA is shown as a color image 
of the W51 IRS2 region. The contours are -14, -13, -12, -11,-10, -9, -8, -7, -6, -5, -4, 4, 5, 6, 7, 8, 10, 11, 12, 13, and 14 times 970 mjy beam"', the rms noise of the image. The integration 
is over a velocity range from 50 to 80 km s"'. The synthesized beam is 3.4 " X 3.2" with a P.A. = -87°, and is shown in bottom left comer . The scale bar indicates the molecular line 
emission and absorption in Jy Beam"' . The white diamonds indicate the positions of the infrared sources KJD3 and IRS2d jKraemer et alJ20QltlOkamoto et alJ2001l ; lLacv et alJ2007[) and 
the radio source W5 1 d2 I Gaume et al. 1 993 ). The yellow line indicates the orientation of the SiO(5-4) high velocity bipolar molecular outflow centered on the water masers (Zapata et al. 
in prep.). The green circle indicates the position of the center of W5 1 North ^Schneps et aU198ii) . 



2. Observations 

The observations were made with the SMAQ during 2005 August 
20. The SMA was in its compact configuration, which includes 
21 independent baselines ranging in projected length from 16 to 
50 m. The phase reference center of the observations was R.A. 
= 19h23m43.80s, decl.= 14°31'30.0" (J2000.0). The frequency 
was centered at 217.1049 GHz in the Lower Sideband (LSB), 
while the Upper Sideband (USB) was centered at 228. 1049 GHz. 

A close blend of the CN N=2-l, 7 = 5/2- 3/2, F = 
5/2 - 3/2 and N=2-l, 7 = 5/2- 3/2, F = 7/2 - 5/2 lines 
were detected in the USB. Their frequencies are 226.874166 and 
226.874745 GHz, respectively. Both lines have very similar in- 
trinsic strengths and energies above the ground state. The LSR 
velocity scale in this paper is given with respect to the rest fre- 
quency of the former line. The velocity difference corresponding 
to the frequency difference is 0.77 km s '. When, in this paper, 
we refer to "the" CN line, we mean this blend. The C2H3CN 
Jk^.k^ - 232,22 - 222,21 line was detected in the LSB at a fre- 
quency of 217.497585 GHz. 

The full bandwidth of the SMA digital coiTelator is 4 GHz 
(2 GHz in each side band). The correlator was configured with 
spectral windows ("chunks") of 104 MHz each, with 128 chan- 
nels distributed over each spectral window, providing a resolu- 
tion of 0.8125 MHz (1.1 km s') per channel. 

The zenith opacity measured {t23ochz) with the NRAO 
tipping radiometer located at the Caltech Submillimeter 
Observatory (close to the SMA) varied during the night between 
0.12 and 0.20, indicating good weather conditions during the 
observations. Phase and amplitude calibrators were the quasars 
1749-1-096 and 1741-038, with measured flux densities and for- 
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mal fitting errors of 2.08 + 0.05 and 1.81 + 0.05 Jy, respec- 
tively. The uncertainty in the flux scale is estimated to be 15- 
20%, based on the SMA monitoring of quasars. Observations of 
Uranus provided the absolute scale for the flux density calibra- 
tion. Further technical de scriptions of the SMA and its calibra- 
tion schemes are found in 'H o et al] (|2004|) . 

The data were calibrated using the IDL superset MIR, 
originally developed for the Owens Valley Radio Observatory 
( IScovilleetal.lll993h and adapted for the SMA0The calibrated 
data were imaged and analyzed in the standard manner using the 
MIRIAD and AIPS packages. We used the ROBUST parame- 
ter of the INVERT task set to -2, which corresponds to uniform 
weighting to achieve the maximum angular resolution while sac- 
rificing some sensitivity. The resulting image rms noise of line 
images was 30 mJy beam"' for each channel at an angular res- 
olution of 3.4" X 3.2" with a RA. = -87°. The data were self- 
calibrated in phase and amplitude using as a model the contin- 
uum image. The final images were shifted ~ 1" in right ascen- 
sion in order to be consistent with the positions of the molec- 
ular cores associated with the W51 North and W51 d2 better 
determ i ned with the Very L arge Array observations of iHo et al] 
(Il983l) : IZhang & Hoi (|1997[) . This discrepancy is mainly caused 
by the baseline error, the finite S/N, and the atmospheric fluctu- 
ations in our millimeter wave observations. 



3. Results and Discussion 

3.1. The Molecular and Millimeter Continuum Emission 

In Figure [1] we present the integrated CN line emission im- 
age of t he W51 IRS2 re gion, and the positions of the H2O 
masers (llmai et al.l |2002|) . the infrared sources W51 IRS2d 



^ The MIR-IDL cookbook b y C. Qi can be found at 
|http://cfa- www.harvard.edu/~cqi/mircook.htmll 
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Fig. 2. Spectrum of the CN N=2-l, 7 = 5/2- 3/2, F = 5/2 - 3/2 and N=2-l, / = 
5/2 - 3/2, F = 7/2 - 5/2 lines towards the center of the W51 North object. The spectral 
velocity resolution is 1.1 km s"^. The dashed line indicates the systemic velocity (VisR=59 
kms"'). 



not allowing us to search for similar east-west velocity gra- 
dients associated with these line molecular tracers. The W51 
No rth sourc e is as socia ted with the extend ed hot core found 
by iHo et all (Il983l) and IZhang et alj d 19981) . and it is centered 
on th e cluster of m asers as reported in o ther obse rvations e.g. 
NH.:lHoetal.l(IT98 3): Eisn er et all (l2002h . CH,CN: IZhang et al.l 
(Il998h . S07: ISollinsetal.l (120041) . If we assume that the molec- 
ular gas is rotating as a rigid body (i.e. the dynamical mass is 
M^y„=v^rsin^(/)/G, where v is the rotation velocity, r is the ra- 
dius of the envelope, / is the inclination angle of the envelope 
assumed to be 90° and G is the gravitational constant), we esti- 
mate a mass for the central object(s) of 40 Mq. This central ob- 
ject might be associated with a single central O-type (proto)star 
or with a small group of B-type (proto)stars. However, as there 
is a strong and compact bipolar outflow in the center of the core 
(see Figure 3), it seems to be dominated by one central massive 
star. 

From Figure 3 and assuming that at a wavelength of 1.3 
mm we are observing isothermal optically thin dust emission 
with a dust mass opacity coefficient that varies with frequency 
as A- oc 1/, with y6 = 1, (the size of the source suggest that 
we are observing emission from the env elope; hence, we adop t 
j0=l, however, this value is uncertain, seel Beckwith et al](ll990h 
that observe how this value varies in pre-main-sequence stars), 
a gas-to-dust ratio of 100 (which may not be the most adequate 
to use for protostellar sources since erosion of the circumstellar 
envelope by photoevapo ration from near OB s t ars may decrease 
the ga s-to-dust ratio, see lWilhams et al.l ( l2005h : lThroop & Ballvl 
(l2005b ). an adopted value of Ki^mm = 1-5 cm^ g~' and a dust 
temperature value of about 100 K (with uncertainties of a factor 
of about 1.5, Z hang et al. (1998) ), we estimate an enclosed mass 
of the molecular cor e W51 North of 90 Mq, very close to the 
100 M0 estimated bv lZhang & Hol(ll997h . Due to the uncertain- 
ties referred to above, the values of the derived masses are good 
within a factor of 2. 



is59kms-' lzhangetal.|[T998l) . the location of the redshifted 3.2. Forming an Early O-type star in W51 Nortli? 



and KJD3/OKYM 1 dKraemer et al.ll2001l; lOkamoto et al. 20011: 
Lacy et al] l2007h . and the radio source W51d2 (IGaume et all 



1993b . all of them in the neighborhood of the W51 North ob- 



ject. Furthermore, we mark the position and orientation of the 
strong high velocity molecular bipolar outflow traced by the SiO 
v= 0; 7 = 5 - 4 line found by Zapata et al. (in prep.), which is 
centered on the cluster of masers. In this image we can see two 
components of the molecular line distribution, one in unresolved 
absorption toward the W5 1 North object and the other also quite 
compact, but resolved in emission surrounding this region. We 
interpret this appearance the signature of molecular absorption 
against the strong compact millimeter continuum source associ- 
ated with W51 North and shown in Figure 3. 

Figure |2] shows the spectrum of the CN line towards 
the center of the W51 North object. The line shows a pro- 
nounced inverse P Cygni profile. The emission feature appears 
at a Vlsr-SS km s"' and absorption at Vlsr=(^5 km s"\ 
which is cons i stent w ith the HCO^ molecular observations of 
lRudolphetal.1 (Il990|). Given th at the core's systemic velocity 



absorption projected against the bright continuum emission of 
the central highly obscured object implies inward motions away 
from the observer. The brightness temperature of the emission 
from the infalling material is lower than ~ 10 K, the bright- 
ness temperature of the cont i nuum emission from the core. From 
Figure 6 of Rudolph et al.l (Il990h and Figure 2 shown here, 
we estimate that the velocity of infall {Vi„faii) is about 4 km 
s ^ With this information and taking the values of the density 
(p=2xl0* cm"^), linear radius (r=1.4 xlO"* AU) rep orted for the 
compact continuum source located in W51 North ( Zhang & Hoi 
[1997; Zhang et al. 1998) and a radius of the hot core of 2 x 10"^ 
AU, and following Beltran et al. (2006), we calculate that the 
mass infall rate (M-Anr^ pV i„faii) is between 4 and 7 x 10"^ M© 
yr"'. This value have large uncertainties, due to the uncertainty 
on the density and on the radius at which V,„/a// is measured. 

The first moment map of the C2H3CN Jk„,k, - 232,22-222,21 
line is shown in Figure 3. The emission is tracing an unre- 
solved east-west rotating molecular "envelope" or "core" with 
a total velocity shift of 1.5 km s"' and a size of 4 x 10'* AU. 
Moreover, the integrated emission from this molecule is well 
centered on the millimeter continuum source, suggesting that 
this species is tracing high density gas close to the (proto)star 
(see Figure 3). Finally, our bandpasses contained lines from 
other molecules (e.g. HCOOCH3 and CH3OH) associated with 
W5 1 North. However, they were very much contaminated by the 
emission from the hot molecular core associated with W51d2, 



The combined 1.3 mm continuum, C2H3CN and CN data from 
W51 North suggest that this object is forming a massive 05- 
type (proto)star in its center through molecular gas accretion and 
with a very large accretion rate between 4 and 7x 10"^ M© yr"' . 
Moreover, the powerful compact bipolar high velocity (> 100 
km s"') outflow traced by the H2O masers pinpoints the position 
of this putative central massive object. 

Zapata et al. (in prep.) in addition found a compact high 
velocity SiO bipolar outflow with both unresolved-lobes spa- 
tially separating for less than one arcsecond and forming a R 
A. of 150° + 30°. This is centered at the positio n of the com- 
pact H2O outflow which has a RA of 1 10-145° dSchneps et al.1 
Il98ll:llmai et al.l2002l;lEisner et al.l2()02h . We propose that these 
two outflows might be manifestations of a single powerful out- 
flow. The masers appear to be tracing the innermost regions 
of it (as observed in other outflows, e.g. IRAS 20126-1-4104, 
Moscadelli et al. 2000), while the SiO is tracing the more ex- 
tended shocked molecular gas. From this point of view, the 
southeastern cluster of H2O masers (see Figure 3) maybe trac- 
ing another older ejected bow-shock. However, higher angular 
resolution SiO molecular observations are neccesary to confirm 
this picture. It is interesting to note that the RA. of the veloc- 
ity gradient across the C2H3CN envelope (~ 90°) is not ex- 
actly perpendicular to the orientation of the molecular outflow 
(RA. ~ 140°), as might be expected. This may indicate that out- 
flow could be precessing due the presence of binary system or 
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At present, there is a list of early massive (proto)stars that 
have been associated with possible infalling motions and with 
large mass accretion rates, e.g. W51e2: a gas mass of 200 
Mq and an accretio n rate of 10"^ Mq yr ' (IZhang & Hdll997l: 
iHo&Yound 119961) : NGC7538-IRS9: a g as mass of 100-30 
Mo and an accretion rate of 10"-^ Mq yr ' JSandell et al.ll2005h : 
G24.78+0.08: a (proto)stellar mass o f ~ 20 Mq and an accretion 
rate between 10"^ to 10 '* Mq yr"' dBeltran et al. 2006): IRAS 
16547-4247: associated with an O -type (proto)star a nd an ac- 
cretion rate of about 10"^ Mq yr ' (" Garav et ad2007') and W51 
North, a (proto)stellar mass of ~ 40 Mq and an accretion rate be- 
tween 4-7 X 10"^ Mq yr"' (these results). This suggests that the 
very massive stars (O-type) might form starting with very high 
accretion rates, sufficient to quench the formation of an UCHII 
region. We note that this hypothesis has already been proposed 
by recent numerical simulations (Ban eriee & Pu dritz 2007). 

We thank the anonymous referee for many valuable sug- 
gestions. R.G. is grateful to the Alexander von Humboldt 
Foundation for a Humboldt Research Fellowship. 
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fig. 3- Upper: SMA 1.3 mm continuum emission color image overlayed with the mo- 
ment zero distribution of the C2H3CN line (pink contours) of the W51 North region. The 
integration is over a velocity range of 55 to 65 km s"' . The synthesized beam with a FWHM 
3.4 " X 3.2" and a P.A. of -87° is shown in the bottom left corner. The contours are -4, 4, 
8, 12, 16, 20, 24, 30, 40, 50, 60, 80, 90, 100, 120, 150, 170 and 200 times 40 mJy beam-' 
km s"' , the rms noise of the image. The scale bar indicates the continuum peak flux den- 
sity in Jy beam"' . Lower: First moment color image of the C2H3CN emission toward the 
W51 North region. The integration is over a velocity range of 55 to 65 km s"'. The scale 
bar indicates the velocity shift in k m s~'. The white diamond s indicate the positions of the 
infrared sources KJD3 and lRS2d fKraemer et al. 2001; Okamoto et al."2001"; 'Lacy et al." 
|2007} and the radio source W51d2 iGaume et al. 1993). The blue and red crosses indicate 
the position of the blue- and red-shifted strong H2O masers spots, respectively, reported by 
[imai et al. 1 2002). Note that the central cluster of masers is tracing a high velocity outflow 
with a northwest-southeast orientation ilmai et al. 2002; Eisner et al. 2002). The yellow line 
indicates the orientation of the SiO(5-4) high velocity bipolar outflow centered on the water 
maser s (Zapata et al. in pre p.). The black circle indicates the position of the center of W5 1 
North ISchnepsetalJ1981h . 



that the C2H3CN emission is contaminated by the outflow. This 
physical phenomenon of precessing outflows has been reported 
in othe r outflows: IRAS 2q i26-H4104 (Shepherd e t al. 200^ 
LI 15 7 (lBachilleretal.ll2001h . and NGC7538IRS1 (iKrausetalJ 
l2006h . 
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